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The assumption of an exact isospin symmetry would imply equal strengths for mirror El transi- 
tions (at least, in the long-wavelength limit). Actually, large violations of this symmetry rule have 
been indicated by a number of experimental results, the last of which is the 67 As - 67 Se doublet in- 
vestigated at GAMMASPHERE. Here, we examine in detail various possible origins of the observed 
asymmetry. The coherent effect of Coulomb-induced mixing with the high-lying Giant Isovector 
Monopole Resonance is proposed as the most probable process to produce a large asymmetry in the 
El transitions, with comparatively small effect on the other properties of the parent and daughter 
levels. 



I. INTRODUCTION 

The presence of symmetries in physical laws in most 
cases greatly enhances our understanding of their nature 
and consequences. Symmetries, either exact or only ap- 
proximate, have a particular importance in the fields of 
elementary particle and nuclear physics. The approxi- 
mate charge independence of nuclear forces, ultimately 
related to the near degeneracy of up and clown quarks 
[l|, permits to treat protons and neutrons as different 
states of the same particle (the nucleon) and to clas- 
sify nuclear states according to the different represen- 
tation of a symmetry group, the Isospin SU(2). In this 
scheme, protons and neutrons are characterised by the 
isospin quantum number T = 1/2, with third component 
T3 = +1/2 and —1/2, respectively. States of nuclei with 
the same mass number A can be grouped, according to 
the value of the isospin T, in isospin multiplets of 2T + 1 
states belonging to the different nuclei, distinguished by 
the value of T3 = (Z — N)/2. Isospin symmetry is vio- 
lated by the electromagnetic interaction (mostly due to 
Coulomb forces among protons) and, to a lesser extent, 
also by nuclear forces. However, the most important part 
of the Coulomb interactions is diagonal with respect to 
T3 and mainly contributes to the mass difference between 
various members of the isospin multiplet. Finer effects 
of the symmetry-breaking forces can be investigated by 
measuring the so-called Mirror Energy Differences [|[ or, 
more generally, differences in excitation energies among 
members of a multiplet. In recent years, this field has 
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become the object of a considerable number of experi- 
mental and theoretical studies, as the level schemes of 
nuclei with T3 = +1/2 (i.e. Z=N+1) could be measured 
for increasingly larger values of A. Furthermore, when 
transition probabilities could be determined, their com- 
parison between mirror nuclei opened an important win- 
dow to investigate the amount and the origin of isospin 
violation. 

Here, we limit our discussion to the relatively simple 
case of El transitions The El transition operator 
is expected to be pure isovector, at least in the limit 
of long wavelengths, where Siegert's theorem Q holds. 
This fact implies that (1) El transitions with AT = 
in nuclei with Z = N are forbidden, and that (2) cor- 
responding El transitions in mirror nuclei have equal 
reduced strength. Both rules arc to some extent vi- 
olated by isospin-non-conserving (mainly, Coulomb) in- 
teractions. In the Z=N case, these violations appear as 
second order effects, while in mirror nuclei the effect is of 
first order. The difference is due to the interference be- 
tween the irregular amplitude (symmetric with respect 
to the exchange of the two nuclei in the doublet) with 
the regular amplitude (which is isovector, antisymmetric 
with respect to the exchange). 

In the following, we discuss the relative importance 
of different possible sources of asymmetry in mirror El 
transitions. As a simple example, we consider in partic- 
ular those nuclei which can be described by the nuclear 
shell model in a limited Hilbert space, containing a full 
major shell and the unique-parity intruder from the next 
major shell. Although the particle-hole excitations in- 
volving all states of the higher shell must be considered 
for a reliable description of the El transitions, we assume 
that the largest part of the El amplitudes only involves 
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the intruder orbital ji and, as a consequence, only the 
largest-j orbital, j'jv = ji — 1, of the lower major shell. It 
is important to note that the inclusion of more orbitals 
in the calculation, briefly discussed in Appendix [Cl does 
not change substantially most of the results. 

Actually, strong asymmetries in B(E1) values have 
been observed in several mirror transitions, e.g. pairs 
of mirror nuclei of the sd and pf major shells The 
clearest examples, however, were found in light N ~ Z 
nuclei, such as 17 (3 and 17 F. Such nuclei often exhibit 
large differences in the neutron and proton binding en- 
ergies, and coupling to the continuum needs to be taken 
into account. The present discussion is limited instead 
to heavier mirror nuclei, in which the smaller binding en- 
ergy of the proton is compensated by the larger coulomb 
barrier. 

As a typical example ("benchmark" in this work), we 
consider the mirror pair 67 As - 67 Se, whose structure in- 
volves the pf shell plus the gg/2 intruder orbital. This 
doublet has been investigated in a recent experiment at 
GAMMASPHERE @. Two pairs of mirror transitions 
with a sizable El component have been observed, con- 
necting the lowest 9/2 + state to lower lying 7/2 _ lev- 
els (Fig. [T]). The measured El strengths and the abso- 
lute value of the corresponding El matrix elements are 
reported in Table [U The 9/2 + state has presumably a 
rather pure gg /2 character, while the daughter states have 
a complex structure and contain only a small component 
that can be reached by the El transition. As a con- 
sequence, the observed values of B(E1) are very small. 
All numerical results reported in the following will re- 
fer to this particular pair of nuclei. The radial integrals 
have been obtained with single-particle wavefunctions in 
a Woods-Saxon potential with spin-orbit interaction, as 
specified in Q. These integrals change slowly with the 
atomic number, and for the / 7 / 2 — > d 5 / 2 transitions in the 
middle of the sd shell they would give results very close 
to those of the g 9 / 2 — > jV/2 transitions in mass A=67. 

In Section|TTl we derive the expression of the El transi- 
tion amplitude from the intruder state a ( J a — ji) to one 
of the normal-parity states b (Jt). No specific assump- 
tions are made on the structure of these states, apart 
from the fact that orbits of the higher major shell, dif- 
ferent from the intruder, give negligible contribution. 

In the following sections, we discuss the different pro- 
cesses that can lead to the presence of an (induced) 
isoscalar El transition amplitude, in addition to the main 
isovector term. In Section IIIII we consider the effect of 
higher order terms in the nucleonic current, in addition 
to those considered in the Siegert theorem, which are in- 
creasingly important when the long-wavelength assump- 
tion fails. 

In Section [TV] we discuss several simple effects related 
to the mixing of wavefunctions: the Coulomb mixing be- 
tween neighboring states (jIV A[) and between states of 
very similar structure, such as the analogue-antianalogue 
mixing PVB[) . None of the processes considered up to 
this point seems able to justify the observed asymmetry. 
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FIG. 1. (adapted from [(|) Partial level scheme of 67 Se and 
67 As, showing the decay of the lowest 9/2 + state. Energy 
labels are in keV. 



We can conclude that the difference in the wavefunctions 
of the two mirror nuclei involves many (weak) mixing 
with a large number of states, possibly lying rather far 
in energy from the levels considered. There are two dif- 
ferent approaches to consider this situation. In the most 
direct treatment , the residual interactions in the two 
mirror nuclei are assumed from the start to be different 
and to include the Coulomb interaction (as well as other 
possible isospin violating terms). It is well known that 
most of the El strength is shifted to higher-lying col- 
lective states, while the low-lying El transitions remain 
substantially hindered, due to the destructive interfer- 
ence among the individual contributions. If the residual 
interactions are not identical in the two mirror nuclei, 
the negative interference can amplify substantially these 
differences in the resulting B(E1). This mechanism is 
easy to understand, but even if a shell model calcula- 
tion in this necessarily huge Hilbert space were to be- 
come possible, the results would scarcely be transparent 
with respect to the nature of the processes involved. One 
could consider, however, the same problem from a differ- 
ent point of view. Namely, let one suppose that a zeroth- 
order calculation were performed with isospin-conserving 
residual interaction. As a next approximation, Coulomb 
interactions could be included to evaluate, to first order, 
the mixing among zeroth-order states. As in the former 
approach, one should expect that coherent contributions 
from collective states play an significant role in producing 
the El asymmetry, as well as the concentration of the El 
strength in the collective states has a role in depleting the 
El strengths of low-lying transitions. The advantage of 
this approach is that it can give semi-quantitative predic- 
tions on the B(E1) asymmetries, even without knowledge 
of their absolute value. Furthermore, it would elucidate 
the principal process (or processes) responsible for the 
largest part of the observed effects. 

A process of such kind, which could in principle ac- 
count for the magnitude of the observed effects (namely, 
the coherent contribution of states belonging to the gi- 
ant isovector monopole resonance) is discussed in detail 
in section HV Cl 
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TABLE I. Values of B(E1) for the transitions proceeding from 
the lowest 9/2 + state in 67 Se and 67 As, as deduced from life- 
times and M2/E1 ratios, determined in [f|. 



Nucleus 


K-.,- 
[keV] 


B(E1) 
[e 2 fm 2 ] 


\(I + \\M(E1)\\1 
[e fm] 


~)l 


67 As 


725 


1.4 ±0.4 x 10 -6 


3.7 ±0.5 x 10" 


-3 


67 Se 


717 


0.4 ±0.4 x 10~ 6 


2.0 ±2.0 x 10" 


-3 


67 As 


319 


8.3 ± 2.5 x 10~ 6 


9.1 ± 1.3 x 10" 


-3 


67 Se 


303 


< 1.4(9) x 10" 6 


< 3.7(11) x 10 


-3 



II. ISOVECTOR AND ISOSCALAR 
CONTRIBUTIONS 

In the following calculations, the El transition is as- 
sumed to take place from an intruder single-particle or- 
bital ji to a normal-party orbital jn = ji — 1 (or vice- 
versa). The parent state a will be the lowest intruder, 
with J a = ji and parity 7f. A possible daughter state b 
must have J b = J a ± 1 or J a and parity tt = —tt. Its 
wave function can contain pairs inside intruder orbitals 
coupled to zero: in this case, the transition could proceed 
from a j'jv orbital present in a to a jj orbital in b. 

If we expand the wave functions of states a and b in 
terms of products of the one-body wave function times 
the core wave function (with the proper fractional parent- 
age coefficients) the only terms of the expansion that con- 
tribute to the transition are those having a common core 
state (of isospin T c = or 1) for both states a and b, and 
the single-particle orbit changing from jj to jm (with 
a core state of positive parity) or vice-versa (with a 
core state J~ of negative parity) : 

\a; J a ,Af a ;l/2,T 3 > = (1) 

J2 C fp {a\jj;^ J+,T C )[<^7) ® *(/*, tf,Tc)]%:% 2) 

± O fp (a\j N ; fx, J-,T c )[cf>(j N ) ® J", T c )}^ 2) 



\b;J b ,M b ;l/2,T 3 > = 

C fp (b\j N ; fx, J+,T c Mj n ) ® *(ft, J+, T c ) 



(2) 



(Jb, 1/2) 
M b ,T 3 



where U^ K > is a tensor operator of rank K acting only 
on the subspace "1" , we obtain for the reduced matrix 
element between states a and b 



(b,J b ,T b \\M^ K) \\a,J a ,T a ) = 

E( 1 \T c +K+lrrt rp /V 2 T b T c 
p [ L) la±b \ T a 1/2 K 

x[Y(-l) Ja+jN+J t +1 J a J b { j j N a f T \ 



x n C fp (a\ji;fx, J+,T C ) C fp (b\j N ; fx, J+ ,T C ) 

.j2(-iy° +ji+j ; +i j a j b / jl Jb J » 



Ja 3N 1 



x n C fp (a\j N ; fx, J ^ , T c ) C fp {b\ji; fx, J^ , T c ) 

x{h\\M% K) \\j N ) 



(4) 



where n is the number of active nucleons, J = y2J ± 1, 
and the triple bars indicate a reduced matrix element 
with respect to ordinary space and to isospin space (with 

T a = T b = 1/2). The operator J\4^ K ^ is now a tensor of 
rank 1 in the ordinary space and K = 1 or in isospin 
space. 
Now, 



{jN\\M^ K) \\3 I )=(-ir-^{ 3 ,\\M^ K) \\j N ) 



= +(ji\\M% K) \\ 3 



N) 



(5) 



Here, ji — ]n = 1, both single-particle states have isospin 
1/2 and the El operator is odd under time reversal. We 
obtain therefore for the reduced matrix element in ordi- 
nary space 



b,J b ;T b ,T 3 \\M ( E K) \\a,J a ;T a ,T 3 
_ r_iy/2-T s I Tb K T, 



-T 3 T 3 



(J 



N 



rib-) 



T C =0,1 k ' 



(0) 



Y c fP (b\j r ,Li, j-^mjj) ® $(^, j-,T c t h ' i/2) 



!M b ,T 3 



Taking into account the relation, to be used both in or- 
dinary space and in isospin space, 

(In ® j2\J\\u (K) {i)\U ® h]J') = (-iy^+ J '+ K 

xV(2J±l)(2J'±l)| j J, j 1 (jiW (K) W > (3) 



where MS 11 } = D^yT and M^ 10 * 1 = are the isovec- 
tor and isoscalar part of the single-particle electric dipolc 
operator: 



(jAHI|A4 10) llb/) = (^P$|j,)(l/2|l|l/2) (7) 



'is i 

UnIM^Ijj) = (j w ||^|b 7 )(l/2||f||l/2) 



and the core-isospin dependent coefficients A(T C ) (T c = 
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or 1) are 

A(T C ) = (-l) Tc 2J a J b 

x \j2(-v J z +Ja+jj 

x nC /p (a|j/;^, J+,T C ) x C fp (b\j N ; p,, J+ ,T C ) 

+ D-^** 4 -* { J /; f } 

x nC fp {a\j N ;fi,J- ,T C ) x C fp (b\ji; fi, J~ ,T C ) 



(8) 



jl 1 



With T a = T b = 1/2, T c = or 1, by inserting the 
numerical values of the coefficient^] one obtains 



(&,J fe ;l/2,T3||^ 11) |a,J ;l/2 ) T 3 ) 

(_l)l/2+T 3 



(9) 



[.4(1) +3.4(0)] (ji\D^\j N ) 



(b 1 J b ;l/2,T 3 \\M i ™ ) \\a,J a ;l/2,T 3 ) 



(10) 



\ [4(1) -4(0)] Ci/PSliiv) 



The leading isovector term in the single-particle operator 
is, in our case, 

(Jill^lliiv) = \ < h\r\j N > (ji\\Y^\\j N ) (11) 

where is the spherical harmonic for £ = 1. 

The different forms of possible isoscalar contributions 
are discussed in the following sections. 



III. HIGHER-ORDER TERMS AFTER SIEGERT 

It is well known that the usual expression of electric 
transition amplitudes, deduced from the Siegert's the- 
orem, is only valid in the long wavelength limit. The 
complete expression for the electric transition amplitude, 
taking also into account relativistic corrections, is given 
by Friar and Fallieros Q: 



,L-1 



T(E,LM) = 



(2L + 1)!! 
x / dr 3 



(12) 



1 Namely, (l/2||l||l/2) = y^; (l/2|r|l/2) = v^; 
( -T 8 T 3 2 ) = (~V 1/2 - T3 /^ ** K = 0, and = 1/VG 

, „ - / 1/2 K To \ 
for K = V > { 1/2 1/2 if J 
= (— l) Tc — 1 /2 for all other cases. 



1/6 for if 



1 and 



where is the vector spherical harmonic and 



9l{z) 
hh{z) 



1 - 



Lz 1 



2(L + 2)(2L + 3) 



(13) 



L + 2 d 

Lz dz 



-2L 



dz 



b 2L+1 9L{z)] 
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For our purposes, it will be sufficient to consider only 
the first term after the Siegert limit, as given in the 
Eqs. (|13[) . We will consider first the part of the integral 
(TT21) containing the time derivative of the charge density 
P = Poif) ex P( — ikct), approximating the nucleus to an 
ensemble of point-like nucleons: 



p = —ikcp(r) 



-ikc^^ j — ^— — — e 5(r — fl 



(14) 



The isovector part of the El transition operator is, for a 
single-particle transition, 

D iv (E1)t 3 = ^|fcc ry«(f)r 3 . (15) 

To first approximation, the isoscalar part of the transi- 
tion amplitude only comes from the second term of the 
series expansion of <7i(fcr). As we have to deal with a 
one-body operator, we can easily obtain the amount of 
this correction with respect to the main (Siegert) term, 
for each single-particle transition: 



(£i,h\Df s (Ei)\£ a ,h) fc 2 (h,h\r 3 \£ 2 ,h) 



30 (h,ji\r\£ 2 ,j 2 ) 



.(16) 



Now we can estimate the numerical value of this ratio 
for the case of the A = 67 doublet [f| chosen as a suit- 
able benchmark, and for aj 9( / 2 -> fj/2 transition. With 
Woods-Saxon radial wavefunctions one obtains 



</7/ 2 |(r/i?o) 3 |39/2) 



= 0.834. 



(f7/2\r/Ro\g 9 /2) 
Assuming R = 1.27-4 1 / 3 fm = 5.158fm 



(17) 



& (/7/2k 3 |ff 9 /2) = l_ {kR ^ 2 (fr/2\(r/Ro) 3 \g 9/ 2) 
30 (/7/2M59/2) 30 (f7/2\r/Ro\g 9 /2) 



1.90 x 10~ 5 (£; 7 [MeV] 



(18) 



The evaluation of the second part of Eq. (TT2"]) can be 
easily performed if we substitute the continuous mag- 
netic density p(f) with that of an ensemble of point-like 
nucleons with spin: 



P(r) = Pn 

i 

1-T3« 



f l + 73(i) 






£% + g P Si 



g n Si >S(f- fi) 



\li + (g P - i - g n )s t r 3 (i)} , 



(19) 
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where fi n = eh/(2M p ) is the nuclear magneton (and M p 
the proton mass). On the basis of Eq. (TTgj) , we can ob- 
serve that the magnitude of this term in comparison to 
the first term of Eq. (TT2"j) is given by 



kec 



Tike 
2M p c 2 



0.53 10 -3 E 1 \MeV\ 



(20) 



Here, we are only interested in the isoscalar part, where 
the contribution of the term s is hindered due to the 
numerical factor g p — 1 + q n w 0.76. The evaluation 



of the matrix elements of s ■ y^' 1 ^ and j ■ 
performed in detail in Appendix |A"1 For our benchmark, 
corresponding to a g Q / 2 — > fr/2 single-particle transition, 
one obtains 



(g9/2ll^l/ 7 /2) _ 2V2 k^ 

(99/2 \\Div I/7/2) 3 ec 



(gg - 1 +J7n) 



V6 



2 x 10" 4 (# 7 [MeV]) 



1 

(21 



if we assume that the above description of the magnetic 
density is approximately correct. 

For 7-ray energies around 1 McV, both correction 
terms are far too small to justify the observed asymme- 
tries in A = 67. 



IV. THE COULOMB MIXING OF WAVE 
FUNCTIONS 

If one takes into account the level mixing due to 
the Coulomb interaction V c , the wavefunction of a pure 
eigenstate |ao > of the charge-invariant Hamiltonian is 
changed into a new one, \a! >. To first order, 



> = \ao > + 



< ak\V c \ao > 
E(ao) - E(a k ) 



\a k > , (22) 



where the sum is extended over all states \a k > having 
the same J™ as | ao > , and which may or may not have the 
same isospin. The El transition matrix element between 
the modified states a' , b' is, again to first order, 



< b'\M{El)\a >=< b \M(El)\a > 



(23) 



E(ao) - E(a k ) 

=< b \M(Ei)\a > + < b \M(El)\a > . 

It was assumed, here, that the Ai(El) operator is pure 
isovector. The ensemble of__the first-order corrections 
(shortly indicated as < ba\M{El)\ao >) transforms as 
an even tensor in isospin space. In the T = l/2orT = 
subspaces, it can be considered as an induced isoscalar 
amplitude. 

If T3 = and the unperturbed states ao, bo have the 
same isospin, the first term of the sum (|23l) vanishes and 



only the induced part contributes. Instead, if T3 = ±1/2, 
the first term is the leading one and the other two are 
only first-order corrections. 

The Coulomb potential can be written as the sum of 
an isoscalar, an isovector and a rank-2 isotensor term: 



Vr 



1 



EE' 



l+T 3 (i)l+T 3 (j) 



1 + g(r(i) • r(j) ) 



(24) 



El^)Eli 



EE e 

1 



2 



[T3(*) T 3(j) 



"(0 • t(j) 



The isoscalar part can be included in the charge- invariant 
Hamiltonian. The matrix elements of the isotensor term 
vanish in the T = 1/2 subspace. They could contribute 
to the mixing with a T = 3/2 state but would produce, 
in any case, equal effects in two mirror nuclei. 

Therefore, any difference between mirror nuclei has to 
be attributed to the mixing induced by the isovector term 



(i) 



1 j^i 



(25) 



where V& is, obviously, a two-body operator. It is pos- 
sible, however, to approximate its matrix elements with 
those of a suitable one-body operator (see Eq. 2-104). 
Actually, the second sum in Eq. (|25]l corresponds to the 
Coulomb potential of a system of A — 1 point-like charges 
e/2 associated to all nucleons j different from the nucleon 
i, and we can approximate it with the electrostatic po- 
tential of an uniformly-charged sphere of radius R, i.e., 
for r < R 



<p c {r) ee Ur/R) 



e(A-l) 3R 2 



2R 



2R 2 



(26) 



(slightly different forms of the function f c will be consid- 
ered in the following). With these approximations, 



VP 



Y,2 T ^ c{n) (27) 

i 

= eT^ c {Q) - |r 3 (0[yc(0) - <p c (r)] 

i 

ee eT 3 ^ c (0) + VP 
and, for the potential tp c of an uniformly charged sphere, 
y(i) = e(A - 1) ^ erf r 3 (i) 



E 



(28) 



The first term of Eq. (|27p is diagonal and does not con- 
tribute to the mixing. The second term is proportional 
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to the isovector monopole operator 

e r?T 3 (i) 



The reduced transition probabilities become, up to first 
order 



(29) 



This result will be exploited again in Section IIV CI 

Actually, the use of a constant charge density inside a 
sphere to evaluate the electrostatic potential ip c is some- 
what inconsistent with the Woods-Saxon distribution of 
matter density assumed to calculate the radial wavefunc- 
tions. Moreover, the tails of these wavefunctions extend 
outside the nuclear radius, in a region where ip c would 
decrease as 1/r. Calculations of the electrostatic po- 
tential for a Woods-Saxon density of charge are given 
in Appendix [Bj For small values of r - i.e., as long as 
the charge density of the Woods-Saxon distribution is 
substantially constant and equal to that of the sphere - 
the values of Vc are equal in the two cases, and the 
differences in the calculated integral are always rather 
small. To obtain the same charge density at the cen- 
tre, the radius R of the uniformly charged sphere must 
take a slightly different value from the parameter i?o of 
the Woods-Saxon distribution. Adopting for the Woods- 
Saxon parameters the values suggested by Bohr and Mot- 
tclson Q, Rq = 1.27A 1 / 3 fm, a = 0.67 fm, for A = 67 
one obtains Rq = 5.158 fm and R = 5.430 fm. 

The matrix elements of are in any case very small. 
To produce a sizable mixing of states, it is necessary that 
the effect be amplified due to some particular conditions. 
This can happen, in particular, (i) when two levels with 
equal J" are very close in energy or (ii) have very similar 
wavefunctions, or (iii) when many different levels con- 
tribute coherently to the mixing. We will consider these 
three cases in the following subsections. 



A. Close-lying states 

The simplest possible case is the mixing of two states 
which lie close in energy. As an example, we can consider 
the El decay of a given state a (of spin J a ) towards two 
states bi, b 2 of equal angular momentum J b , and rather 
close in energy. In this case, taking into account only the 
Coulomb mixing between b\ and b 2 (and neglecting small 
isoscalar terms in the El operator) we obtain up to first 
order 

{b' 1 ,J b \M{El)\a,J a ) = (30) 
= (&!, J^M(El)\a, J a ) + aT 3 (b 2 , J b \\M(El)\\a, J a ) 
(b' 2 ,J b \\M(El)\\a,J a ) = (31) 
= (62, J b \\M(El)\\a, J a ) - oO , 3{b t ,Jb\M{E\)\a, J a ) 

with 

aT 3 = < h\VP\b 2 > I [E(bi) - E(b 2 )} (32) 

In fact, as a consequence of the Wigner-Eckart theorem, 
the matrix element of Vc must be proportional to that 
of T 3 . 



B(El;a-> b' x ) 



1 



(b u J'\M(El)\a,J) 2 (33) 



2J+ 1 

-2aT 3 (b u J b \\M{El)\\a, J a )(b 2 , J b \\M(El)\\a, J a ) 
1 



B(M;a-> b' 2 ) = 



{b 2 ,J b \\M(El)\\a,J a y (34) 



2J a + 1 

-2aT s (b u J b \\M{El)\\a, J a )(h, J^M{El)\a, J a ) 
Hence, the sum of the two reduced strengths, 
B(£l;a-> b[) +B(M;a-> b' 2 ) = (35) 
— ^— [(&!, J b \\M(El)\\a, J a ) 2 + (62, J b \M(El)\\a, J a ) 2 

a + 1 - 

is independent of T 3 and consequently identical in the 
two mirror nuclei. If one of the two unperturbed tran- 
sition strengths (either for a — > b\ or a — > b 2 ) is much 
smaller than the other, a large percentage difference be- 
tween mirror values can be found, but only for the weaker 
transition. 



B. Analogue — antianalogue mixing 

A second interesting case concerns the mixing between 
two very similar wavefunctions, as for a pair of analogue 
- antianalogue states (this would be a very favourable 
case of the mixing of T = 1/2 and T = 3/2 states dis- 
cussed in Q). Let us consider, as a simple example, the 
state obtained with the coupling of a jj =9/2 nucleon 
to the lowest state </>o (J 71 " = + , T = 1) of the isospin 
triplet A = 66. Isospin 3/2 states are obtained in the 
two \T 3 \ =3/2 nuclei. In the \T 3 \ = 1/2 nuclei 67 As, 67 Se 
two independent wavefunctions will result from the cou- 
pling, and two pure isospin states can be constructed by 
proper linear combinations: a T = 3/2 state I03 >, which 
is the isospin analogue of those in the IX3 1 = 3/2 nuclei, 
and a T = 1/2 state |ai >, sometime referred to as the 
anti-analogue of them. Here we will give the results for 
the T 3 = +1/2 nucleus (from which, those for T 3 = —1/2 
can be easily deduced by means of the Wigner-Eckart 
theorem): 

|oa) = \[<f>j(t = 1/2) <g> $ (T C = l)]j,T = 3/2) (36) 
= ci\^(g 9/2 ) $ (T 3 =0)) 
+ c 2 \(j) v {g 9/2 ) $ (T 3 = 1)) 

|oi) = = 1/2) 8 MTc = l)]j,T = 1/2) (37) 
= c 2 1^(39/2) *o(T 3 =0)) 

-Cl|<^(<?9/ 2 ) $0(23 = 1)) 

where, for T 3 = +1/2, 

d = (1/2, 1/2, 1,0 I 3/2,1/2) 

= -(1/2,-1/2,1,1 I 1/2,1/2) = V2/3 (38) 
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and 

Pa = (1/2, -1/2, 1,1 | 3/2,1/2) 

= (1/2,1/2,1,0| 1/2,1/2) = VV3- (39) 

We now use Eqs. (|27I28|) to approximate the non diagonal 
part of the isovector Coulomb interaction V^P with a 

one-body operator V@ , whose matrix element between 
analogue and antianalogue states is, for T3 = +1/2, 

= cic 2 (^$ (l,0)|V^ 1) |^*o(l,0)^ 

■c 2 c 1 (^$o(l,l)|^ 1) |^*o(l,l)) 

= ClC 2 



\% l \^)-U\% 1 \^)\<j> 



- (* (l,l)|vW|$o(l,l) 



(40) 



The diagonal matrix element of the isovector operator 
Vq over the core state T = 1, T3 = is zero. 

Starting from Eq. (|40[) and assuming an energy spacing 
E{a^) — E(a\) = AE, we can now estimate at least the 
order of magnitude of the mixing coefficient. For T3 = 
+1/2, 



V, 



(1) 



c 



a 1 



(-AE) 
$o(l,l 



cic 2 - 



8R AE 



2^59/2 1 -^159/2 



$0(1,1; 



(41) 



In the second term, the contributions of a proton and 
of a neutron in the same orbit cancel one another, due 
to the opposite eigenvalue of T3. There are, however, 
two excess protons in the T3 = 1 core state. If all the 
radial wavefunctions of active nucleons in the core were 
equivalent to that of the ji orbit, the second term in the 
sum of Eq. (|4"Tj) would exactly cancel the first one. We 
can expect, therefore, a resulting matrix element sub- 
stantially smaller than the first term alone, due to the 
effect of the core term. However, the expectation value 
of r 2 1 R 2 in the ji = gg/2 orbit is certainly larger than 
those for the lower orbits in the core. For A = 67, and 
with Woods-Saxon wavefunctions, the radial integral of 
(r/R) 2 in the Ogg/2 orbit is 0.7495, while in the normal- 
party orbits 0/7/2- 0/ '5/2, IP3/2 and lpi/2 is, respectively, 
0.6251, 0.5922, 0.6251, and 0.6359. In Eq. gl]), we will 
use the average of these values, < r 2 /R 2 >= 0.6119, and 
the above estimate of the matrix element in the (?g/ 2 or- 
bit, to evaluate an order of magnitude for the analogue- 
antianalogue mixing. Numerically, with cic 2 = V%/3, 
A = 67, R = 4.43fm and assuming AE m 4 MeV as in 
we obtain a w 0.071. As the matrix element of 



2 See Table [TT] in the Appendix [B] With a Woods-Saxon charge 
distribution, the estimate does not change more than a few per- 
cent 



the isovector interaction V(p between a state of isospin 
3/2 and a state of isospin 1/2 is 

(3/2,T3|^|l/2,T 3 ) = (-l) 3 / 2 -^( 3 g I V 2 

x (3/2||V^ 1) ||l/2) (42) 

the value of a has equal sign in both nuclei of the doublet. 

The El transition matrix element from the state \a[ > 
to a given state \b > will be, at the first order, 

(b\\M(El)\\a[) = (6| J M(M)|a 1 )+a(6| J M(M)|a 8 >.(43) 

We assume, for sake of simplicity, that the state b has 
pure isospin 1/2. If, as we have supposed, the El tran- 
sition proceeds from a j j = 9/2 to a jjy = 7/2 single- 
particle state, we can use for the state \b > a fractional 
parentage expansion in the style of the first line of Eq. © . 
But only the terms corresponding to the coupling of a 
nucleon in the state jm = 7/2 to the core configuration 
$0 with J = 0, T = 1 can be reached by the El transi- 
tion. We can write the (presumably small) part of the 
wavefunction of the state b which is relevant for the El 
transition in the form of Eq. (|3"T|) . 

\b) = |[&(t = 1/2) ® $o(T c = l)]j,T = 1/2) + ... 
= c 2 |<M/ 7 /2) $o(T 3 = 0)) (44) 
-cilM/7/2) $o(T 3 = 1)> + ... 

Taking into account the effective charges for El transi- 
tion, e n = 1/2 and e„ = —1/2, from Eq. (|4"3")l we obtain 



(6|.M(£;i)|K) = [(|c 2 | 2 e w + | Cl | 2 e,) + a Cl c 2 (e n - e„)] 
x (f 7 /2\\erY^\\g 9/2 ) 
-1 + 2V2 a. 
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(fr/2\\erY^\ 



99/2) 



For T 3 = ±1/2, using Eq. ([42j) we obtain the numerical 
coefficient (+l+2-\/2 a)/6. In conclusion, the El strength 
in the two mirror transitions is proportional to (+1 + 
2\/2a) 2 . The mirror asymmetry in the El strength is 
therefore, approximately, 



B(E1, As) - B(E1, Se) 4^2 



a 



B(El,As) + B(El,Se) l + 8a 2 



0.386 



(45) 



and the ratio B(E1, As)/B(El, Se) w 2.26. We note, 
however, that such a large asymmetry has been obtained 
for a pure configuration of the analogue and antianalogue 
states, while the antianalogue strength is usually spread 
over a number of final states [llj | , a situation which will 
strongly reduce the mirror asymmetry in the El strength. 
A detailed shell-model investigation would possibly eluci- 
date the role of the analogue-antianalogue mixing in the 
El asymmetry between mirror nuclei, as the analogue 
and the antianalogue states can be described in the same 
shell-model space. 
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C. Coherent enhancement of induced isoscalar El 

The Coulomb mixing discussed in the previous subsec- 
tions involves states belonging to the same set of shell- 
model orbits necessary for the (unperturbed) parent and 
daughter state of the El transitions (presumably limited 
to two major shells). However, it is well known that a 
comparatively large contribution to the isospin mixing 
comes from states outside this model space, as those be- 
longing to the giant isovector monopole resonance [l2j]. 
Obviously, the mixing with any of these higher-lying 
states, induced by the isovector part of the Coulomb in- 
teraction, is expected to be very small. The combined 
effect of many higher-lying states on the El transition 
amplitude can however become appreciable if their in- 
dividual contributions combine coherently. We shall sec 
how this can be the case. 

We have seen (Eqs. I27I28|) that the non-diagonal 
isovector part of the Coulomb interaction V c can be ap- 
proximated with a one-body operator Vq 1 , having the 
same form of the isovector monopole operator (EO). 
Therefore, it is a sensible approximation [l3|, [l4| to con- 
sider in the ensemble of states a k , b k (with k ^ 0) of 
Eqs. (|22I23[) only those of the isovector monopole reso- 
nances built over ao and 6o, and to use the mean exci- 
tation energy AE a (or AE^) of the giant resonance over 
the state ao (or &o) m the place of those of individual 
states. In this case, Eq. ([2"3]) becomes 



< b'\M(El)\a' >= (46) 
=< b \M(El)\a > + < b Q \M(El)\a > , 



where 



<bo\M(El)\a Q > (47) 
]T < b \M(El)\a k >< a k \V c \a > 



AE a 
-1 



< bo\V c \b k >< b k \M(El)\a > 



We are only interested in the isoscalar part of M(El), 
which results from the isovector part Vc of the Coulomb 
interaction. Approximating the non diagonal part of Vc 1 ^ 
with the one-body potential of Eq. (|28[) , the closure ap- 
proximation gives 

J2 < bo\V L (1] \b k >< b k \M{El)\a > 

« < b \vj 1} M{El)\a > (48) 
Y / <bo\M(El)\a k ><a k \vP\ 

m < b \M{El) K (1) |ao > ■ 
and therefore (as A4(El) and Vc commute) 

< 6ol #»(Ei) K >» 2 «a|^ m 



where we have assumed AE a ps AE b 

2 



AE . 



<6 |^ (0) (^l)K (1) |ao >f 



5 "| E |n^ (1) (ri)r 3 (i) £ ItVcfo) - ^ c (0)]r 3 (j) 



(-AE ) 

5> 



a 



bo\M { °l b + M { °l b \a 



where M-fl b is the one-body operator resulting from the 

term with j = i in the second sum, and M-n-b * s a two- 
body operator resulting from all other terms. As r| = 1, 
the first term is 



(50) 



(0) " a^-E elvcW-Vcin)} fnY^Ksi) 



With the expression of (p c corresponding to the uniformly 
charged sphere, given in Eq. (|26p (and extrapolated also 
for r > R) one obtains for the one-body operator 



R 2 2 



(52) 



which has the same structure as the one coming from the 
second-order term in the series expansion of Eq. (|12l) . 
with a different coefficient, 



C = +{A- l)e 2 / (8RAE 



(53) 



An alternative calculation using a Woods-Saxon charge 
distribution is reported in the Appendix [Bl 

As the one-body operator (|49|) is isoscalar, its ma- 
trix elements can be expressed in the form anticipated 
in Eq. ([ID]): 



M4> T 3|El|^y«(n) 



i[A{l)-A{Q)] e - (j W |ryW(f)|iz) 
„ (j N \(r/R) 3 \ji) 



<jN\(r/R)\ji) 



(54) 



Again, we can evaluate the numerical results for our 
benchmark doublet. For A = 67, we assume R = 5.30 fm. 
The energy difference is AEq « 20 MeV in 60 Ni (accord- 



ing to [15|). As AEq is expected to scale as A 1/3 [12j, 
we assume AE rj 19.3 MeV for A = 67. With these as- 
sumptions, the numerical value of the adimensional coef- 
ficient C in Eq. ([5U]) is C = 0.116. For the ratio of radial 
integrals (last factor of Eq. (|54p). with the radial wave- 
functions corresponding to the Woods-Saxon potential 
one obtains (g 9/2 |(r/i?) 3 |/ 7/2 ) / (g 9/2 \r/R\f 7/2 ) = 0.752. 

It remains to consider the two-body term (second term 
of Eq. ([50]) ). Again, we can use the fractional parentage 
expansion of Eqs. ([US). Here, however, the tensor oper- 
ator is the product of two factors: a vector isovector term 
acting on the single-particle state and a scalar isovector 
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one acting on the core state. The product T 3 (i)T 3 (j) con- 
tains an isoscalar and an isotensor part: 



r 3 (i)T3(j) = Mi)T 3 (j)- - (f(i) • r(i)) 
+ 3 • r(i)) 



(55) 



but only the isoscalar is effective if the states clq and bo 
have T = 1/2. To evaluate the reduced matrix element 
for the isoscalar part of the two-body operator 



M 



(0) 
2-b 



6A£ 



(56) 



we can use the standard relations of tensor algebra for the 
matrix elements of tensor products to obtain the reduced 
matrix element (in ordinary space J5 



(bo; Jb\ 1/2, 



T, 



M 



(o) 

2-b 



a ; J , 1/2, T s 



U'iv|^ (1) (r)|j/) 



7! 



^ Tc+1 f 1/2 T c 1/2 

T C ,T' 



(57) 



1 JaJb / Jft jjV Jj~ 



C/p(o|ii; M) J p ) T c) /A J^, t' c ) 

x (M,J+,r c |^^f( J )| Ai ',j;,T ( ;) 



J" I JJV Ja 1 



Ec- 1 )' 1 " 

Cf P (a\j N ; (i, J' , T c ) C fp {b\ji; fi' , J~ \T' C ) 



X ifi, 



as J M ' = J M . As the A^-b operator transforms as a scalar 
in isospin space, its matrix elements have the same sign 
in both nuclei of the isospin doublet. 

The parent state can have T = or 1, and in princi- 
ple we have to consider both diagonal and non-diagonal 
matrix elements (in the parent-state variables) of the 
isovector operator ~}2(r 2 /R 2 ) f(j). Obviously, its matrix 



3 In fact: (1/2, T 3 ; 1/2, T 3 |(r(i) • r(j)\l/2, T' c - 1/2, T 3 ) 
(-l) Tc+1 {^ y 2 X { 2 } (l/2||f(i)||l/2)(T c ||T(i)ll^); and 
Un, J" M , J b \\nY l ~ 1 \(i){r 3 /R) 2 )\\j I .J^ J a ) = J a J b 

31 J» Ja ){j N \\r l YW{i)\\j I ) (JMrj/Rn^) 
10 1 



elements vanish when T or T" is equal to zero. Other- 
wise, we can use again a fractional parentage expansion. 
Only terms having the same parent can contribute to 
the matrix element and, in addition, the one-body oper- 
ator has non-diagonal terms only between single-particle 
states (with equal f) differing by at least two units of 
the principal quantum number: i.e., it does not possess 
non-diagonal matrix elements inside our model space. As 
for the diagonal ones, shells (or sub-shells) completely 
filled with protons and neutrons do not contribute to the 
sum, as they necessarily have T = 0. If the valence nu- 
cleons are all in the same subshell (or, approximately, 
in subshells with similar < (rj/R) 2 >«< {r/R) 2 >„), 
the integral over the radial coordinates can be factorised, 
^ r(i) = 2T, and only the diagonal terms with fx' = fx 
survive. Therefore, the matrix element takes the form 



(n,J„T c j\^(r 2 /R 2 )T(jW,J,X) 



(58) 



« < (r/R) 2 > v (J„|1|J M ) (T C \\2T\\T' C ) 8^, 
where < (r/R) 2 > v is the average over active valence 



nuclcons, and (J/Jl]! JmJ 



Ju 



For T r 



TL 



(T C \\T\\T C ) = y/T c (T c + 1)(2T C + 1) = \/6. By compar- 
ing the result with Eq. (|8]), we obtain approximately (as 
the first 6- J coefficient has the value —1/3): 



b Q ;J b ;l/2,T 3 



M 



(o) 

2-b 



a ; J a , 1/2, T a 



« -C 1 A(l) e -( ]N \rY^(f)\ 3 i) ((r/Rf) v (59) 

Actually, the expectation values of r 2 /R 2 for the differ- 
ent orbitals of the pf shell (estimated with Woods-Saxon 
wavefunctions) do not differ more than 3% from their av- 
erage value 0.615, as we obtain in Appendix [B] By using 
this average value, one obtains for the numerical coeffi- 
cient of the 2-body term (2/3) ((r/R) 2 ) v = 0.410. As this 
value is not negligible in comparison to that of the 1-body 
term (0.752), a sizable quenching of the isoscalar transi- 
tion amplitude corresponding to the 1-body term results 
from the negative interference of the 2-body term. A sim- 
ilar effect is found for the El transitions with AT = 
in the N = Z nuclei fl3| . However, in the present case 
the quenching only concerns the parent T = 1 term. As 
the parent T = term of Eq. has no counterpart 
in the 2-body matrix element, its contribution remains 
unaltered. 

If we assume that the most important contribution to 
the asymmetry is due to the effect of coherent mixing, as 
approximated in this paragraph, we obtain 



<T 3 ) 



(b, 7/2-; 1/2, T 3 \M% 0) R 9/2+; 1/2, T 3 ) 
(b, 7/2"; 1/2, TsiM^ja, 9/2+; 1/2, T 3 ) 



(60) 



( _ 1) i/2+t 3 3C 



(jN\(r/R) 3 \ji) nA(l) - A(0) 
(.]N\r/R\ji) .4(1) +3.4(0) 

where the quenching factor 7/ takes into account the neg- 
ative interference with the two-body term of Eq. (1501) . 
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Equation (|60l) only gives an approximate estimate of 
the effect, due to the many simplifying assumptions (no- 
tably, the closure approximation) that have been intro- 
duced to obtain this result. Moreover, inclusion in the 
model space of other orbitals of the upper major shell 
(as discussed in Appendix [C]) would somewhat alter this 
result. However, it could be instructive to evaluate some 
numerical results, also in the limited space considered, 
to show that the coherent mixing with the IVGMR can 
explain the large values of the El asymmetries observe 
in our example of the A = 67 doublet, while the simplest 
processes discussed in the previous sections were not able 
to do. 

With the above estimate, r) = (0.752 - 0.410)/0.752 = 
0.458, and the asymmetry ratio for the mirror El 
strengths is 



n ^ B(E1,T 3 = -1/2) 



B(E1,T 3 



-1/2) 



1 



1 



(61) 



where we have put e~ = e(— 1/2) = — e(+l/2). Now, 
to obtain a more accurate estimate one should know the 
ratio 4(0)/4(l), which in turn depends on the Cf p coef- 
ficients. 

The relative sign of .4(0) and 4(1) depends on the com- 
bined effect of all terms in the sum of Eq. (|8]). However, 
we can notice that each of them contains a factor (— 1) T . 
If any of these terms dominates, the relative sign of 4(0) 
and -4(1) is well defined and negative. Actually, this is 
very probably the case also under somewhat broader con- 
ditions. Most probably, the second line of Eq. flSJ) (cor- 
responding to negative-parity parents) is only a small 
correction in comparison to the first one. Let us con- 
sider, from now on, the numerical values corresponding 
to the A = 67 doublet. We can note that the expression 



i-iy 



9/2 9/2 J 
7/2 7/2 1 



has always the same (negative) sign for all J values (from 
to 7) and its value changes very slowly as long as J < 3. 
Therefore, unless the parentage coefficients have a very 
singular behaviour, the relative sign is determined only 
by the factor (-l)J (see also Eq. (g))). 

To obtain just an order-of-magnitude estimate of the 
expected effect, we could evaluate the asymmetry in the 
A = 67 doublet, for two limiting cases in which one 
of the two coefficients 4(1) and 4(0) is negligible in 
comparison to the other. Neglecting 4(1) one obtains 
er ps — 0.753C ps —0.0872 and the asymmetry ratio 
K ps 0.705. 

Taking into account also 4(1) would bring to smaller 
asymmetry (larger 1Z) if 4(1) and 4(0) have the same 
sign, but can also result in a larger asymmetry if - as it is 
most probable - they have opposite sign. If, instead, 4(0) 
is negligible in comparison to 4(1), e~ is positive and 
its value depends on the coefficient n, which takes into 
account the negative interference of the core terms. With 
i] = 0.458, for 4o <C 4i one obtains e~ ps +0.120 and 



7Z ps 1.62. Again, a larger asymmetry could be obtained 
if also a contribution from 4o (having opposite sign) is 
included. 

These results do not change appreciably if one assumes 
a charge distribution of Woods-Saxon shape (Appendix 
E]): one obtains n = 0.445; for 4i < 4 e~ « -0.0852 
and Tl ss 0.710; for 4i > 4 , e~ « +0.116 and TZ ss= 
1.58. 

A last comment concerns the expected sign of e~ . If 
the dominant term in the Eq. (|60[) is the one with T = 1 
parent, e~ > and the reduced strength should be larger 
in the nucleus with N = Z + 1 , for all transitions between 
<?9/2 and fj/2- The opposite is true if the T = parent 
dominates. Again, qualitative considerations can help in 
predicting the relative importance of the two terms. It 
is likely, in fact, that one of the most important parents 
be the lowest J — 0. Now, if A = An + 1, the lowest 
J = parent state is the ground state of the even even 
self-conjugate nucleus with A — 1 nuclcons. Instead, if 
A = An — 1 (as in the case 67 As - 67 Se), the selfconjugate 
parent nucleus is odd-odd and the lowest J = parent 
has T = 1. If this consideration is correct, the predicted 
sign of the asymmetry is consistent with the experimental 
results in the A = 67 mirror pair. 



V. CONCLUSIONS 

It seems worth summarising the results obtained for 
the different processes which could, in principle, produce 
an asymmetry in the El transition strength, as observed 
in the case of the 67 As — 67 Se mirror pair. Higher-order 
terms, either of "electric" or "magnetic" origin, usually 
excluded from calculations by the approximation linked 
to the Siegert's theorem, in the case considered are three 
orders of magnitude lower than the leading one. We note 
that these corrections apply to the transition operator 
and not to the level wavefunctions. Therefore, as long 
as - as it was assumed here - most of the shell-model 
terms contributing to the El transition involve the same 
pair of single-particle states, the same combination of 
fractional parentage coefficients is involved for both the 
isoscalar and the isovector term. Thus if the isovector 
term is hindered as a consequence of accidental cancella- 
tion, a similar hindrance factor can be expected also to 
the isoscalar, leaving the ratio almost unchanged. Only 
meson currents, neglected in our approximate estimation 
of the magnetic term, could break, to some extent, the 
above conclusion. 

The Coulomb interaction, mixing in a different way the 
level wavefunctions in the two mirror nuclei, is presum- 
ably at the origin of the observed asymmetries. Its effect 
could be enhanced when a pair of levels having equal J" 
lie, accidentally, close together. E.g., this could have 
been the case for the two 7/2~ levels lying between 640 
and 1100 keV in 67 As and 67 Se. However, if the asymme- 
try originated uniquely from the mixing between the two 
daughter levels, the total sum of the reduced strengths of 
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the El transitions feeding these levels ought to be equal 
in the mirror nuclei, in contrast with the experimental 
evidence. 

The Coulomb mixing could also be enhanced if it took 
place between states with two "very similar" wavefunc- 
tions. In Section IIVBI we considered an hypothetical 
mixing between a "isospin analogue" state and its cor- 
responding "antianalogue" . In the case of mass A = 67, 
this mixing would lead to an asymmetry similar in size to 
the observed effect. It would also give the right sign for 
the asymmetries. However, this would only happen if our 
T = 1/2, = 9/2+ state would be the exact antiana- 
logue of the lowest T = 3/2 state with the same J 17 , while 
some spread of the antianalogue strength among different 
levels is expected also in this region of nuclei [ljl [I?} ■ 

The effects of Coulomb mixing considered thus far only 
involved states in the same Hilbert subspace needed to 
describe the parent and daughter states of the El transi- 
tion: in the simplest full major shell and at least 
one particle-hole excitation to the next major shell. A 
shell-model calculation in this Hilbert space could treat, 
on the same footing, both the regular (isovector) part of 
the El transition amplitude and the "induccd-isoscalar" 
term originating from the mixing. In such a calcula- 
tion, the isovector part of the two-body Coulomb inter- 
action could be added to the empirical residual interac- 
tions, which could also include the symmetry-violating 
part necessary to account for the Coulomb Energy Dif- 
ferences [l8| . 

Finally, we have considered the possible effect of mix- 
ing with states outside the truncated shell-model space, 
as those belonging to the Giant Isovector Monopole Res- 
onances. With the approximations discussed in Section 
IIV C[ this effect could also be expressed in a form that 
could be treated in the truncated space, if the mean ex- 
citation energy of the monopole resonance were at least 
approximately known. 

A shell-model calculation in such a restricted basis 
could therefore be able to identify the origin of the ob- 
served asymmetry in El transition strengths. At the mo- 
ment, the coherent contribution of states belonging to the 
Giant Isovector Monopole resonance appears as the most 
probable candidate. 
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Appendix A: Evaluation of the reduced matrix 
elements for the magnetic term 

Here we evaluate the reduced matrix elements of the 
operators entering in the second line of Eq. (|12[) . between 
single particle states iuji and £ 2 ,j 2 . To this purpose, 
the following property [19| of vector spherical harmonics 
is exploited: 



y 



M 



M 



(Al) 



where v is a generic vector. Here the cases v = s and 
v = j are considered. 

In the first case, the reduced matrix element of the ten- 
sor product can be obtained easily, because s and Y"( L ) 
operate on different Hilbert spaces 



4, Jill 



y^W 1 ) 



2, J"a 



(A2) 



r h 1/2 Ji) . w . 

= JiJJ 2 < h 1/2 J 2 W*i|y (i) |* 2 ) (l/2||s«||l/2j 

where J = L and ( 1 / 2 1| s < 1 > || 1 / 2) = ^/ij2. The relation 
(ii,Ji\Y^U 3 ,J 3 ) =(-1)^+^+^+1/2 (A3) 

can be exploited to express the result in function of 
(4, Ji\\Y^\\£ 2 , J 2 ) as in Eq. dHJ: 



h,Ji 



y(£) 



(L) 



2,J; 



(A4) 



(_l)-fc+'i+W2 y/2 L 



Ji J 2 L 



h 1/2 Ji 
l 2 1/2 J 2 



1 L 



e 2 h 1/2 



h,J\\Y^\l^J^ 



The second case is not so simple, because the operators 
j and y( L ) do not commute, so that the symmetrised 
form of the operator must be employed. Furthermore, 
they operate on the same Hilbert space, but one can ex- 
ploit the fact that j has no matrix elements between dif- 
ferent single-particle states to obtain: 



\ 1 



(■/) 



jW^yW \\£ 2 ,J 2 



Uii,Ji\YW\i 2 ,j2) (-iy h+J+J u 



where (j||j«||j) - + l)(2j + 1). 



(A5) 
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In the present case, L = J = 1, l\ = 4, J\ ~ 9/2, 
^2 = 3, J 2 = 7/2. With these numerical values, the 
coefficients of the reduced matrix e leme nt of y W in the 
Eqs. ([M1E5|) . are, respectively, yA~/6 and -y/T/2. 



Appendix B: Radial wavefunctions and Coulomb 
potential with a Woods-Saxon distribution 

The radial wavefunctions have been calculated assum- 
ing a Woods-Saxon potential plus spin-orbit: 



V(r) 



v + v s e-s^ — 



d 

7h- 



i 



1 -f e (r-Ro)/a 



(Bl) 



with the values of the constants consistent with Bohr 
and Mottelson ||: V Q = -51 McV, V s = 22 McV, R = 
r A (1 / 3) , r = 1.27 fin and a = 0.67 fin. 

For a consistent evaluation of Coulomb interactions, 
one needs the average electrostatic potential (p c (r) of a 
distribution of A — 1 point charges e/2, which will be ap- 
proximated with a continuous charge distribution having 
a Woods-Saxon shape: 



p e (r) 



Po 



where 



Po = (A-lf- 



1 



With the condition that tp c (r) 
obtain 



4nr dr 



for 



(B2) 



(B3) 



co, we 



Po 



Po 



dy 

V 2 Jo l + e" 



x 2 dx 



(B4) 



This integral has been evaluated numerically, for A = 67, 
with the parameter values suggested in 0. In figHJ the 
result is compared with the potential of an uniformly 
charged sphere of charge density equal to po and total 
charge (A — l)e/2. The radius R of the sphere is deter- 
mined by the condition 



47T - 



2po 



1 



47rr 2 dr (B5) 



To simplify the comparison of the results, ip c (j) is ex- 
pressed in terms of the adimensional function f c (r/R): 



AR 



Mr/R) 



(B6) 



and we define A/ C (r) = / c (0) — f c {r). For the (extrap- 
olated) potential of the uniformly-charged sphere, one 
obtains A/ C (r) = (r/R) 2 . One must now calculate the 
matrix elements of the operators and M-^-b 

fined in Section llV CI For the one-body term, we consider 
the ratio 



(f7/2\(r/R)f c (r/R)\g 9/2 ) 
(f7/2\r/R\g 9 /2) 



(B7) 



10 



<Pc(r) 

[e fm" 1 ] 





utj{r) 

[fm- 3 / 2 ] 

0.5 



-0.5 





a 


1 









r [fm 



10 





b 











r [fm 



10 



FIG. 2. (color online) a: electrostatic potential for uniformly 
charged sphere (dashed) and with a Woods-Saxon charge dis- 
tribution for A — 67 (continuous line) . The dotted line shows 
the continuation outside the sphere of the expression for the 
uniform distribution (dashed curve) in the internal region. 
b: examples of radial wavefunctions for the Woods-Saxon po- 
tential (+ spin-orbit) with the parameters suggested in 
099/2 (continuous line), Of 5/2 (dashed), lpi/2 (dotted). The 
vertical dotted line corresponds to the value of the nuclear 
radius R. 



while for the two-body term (and also for the calculations 
of Section liVBp . it is sufficient to evaluate the diagonal 
matrix elements of f c (r/R) 

By numerical integration, with the parameters of Q 
one obtains the values of the necessary integrals reported 
in the last column of Table [TTJ In the other columns, 
the corresponding values are calculated, with the Woods- 
Saxon wavefunctions, for the potential of the uniformly 
charged sphere and for the extrapolation of the inner 
potential outside the sphere (dotted line in Fig. Eh). 



Appendix C: Effect of the inclusion of more orbitals 

Until now, we have assumed that only the intruder 
orbit c/9 /2 is significant for the description of the relevant 
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TABLE II. Values of radial integrals for different assumptions 
on the Coulomb potential. In all cases, < f7/2\r/ R\g$/2 >= 
0.8285. 

Constant p c Woods-Saxon 
sphere extrapol. distribution 



</7/2k A /c|99/2 


> 


0.700 


0.752 


0.739 


</7/2M99/2> 












< 39/2|A/ c |3 9/ 2 


> 


0.697 


0.749 


0.735 


< fr/2\Af c \f r/2 


> 


0.594 


0.625 


0.620 


</5/a|A/c|/ 6/a 


> 


0.564 


0.592 


0.587 


< P3/2|A/ c |p 3 /2 


> 


0.572 


0.625 


0.608 


< Pi/a|A/ c |pi/ a 


> 


0.580 


0.636 


0.617 



states. As a consequence, only the transitions between 
gg/2 and f-j/2 contribute to El. If other orbitals of the 
upper major shell {e.g. 1^5/2) arc taken into account, 
other orbitals of the lower major shell can be involved 
in the El transitions. We consider now the changes that 
must be introduced in our calculations as a consequence 
of the inclusion in the model space of the two complete 
major shells. 

Equation (JTJ must be modified as follows: 

\a;J a ,M a ;l/2,T 3 > = (CI) 

E E C fp(°\3i\ M, Jt,T e )mi) ® 4, T c )]\ J - 1/2) 



iM„.T 3 



Ji M 



and similarly Eq. ([2]). Equation ((BJ becomes 

(6, J b ;T b ,T 3 \M% K) \a, J a ;T a ,T 3 ) = (-l)V»-*i 

T fc A' T a 
-T 3 T 3 



1 Jn 

i) K E A,^m.){!/ a 2 ^ J} (C2) 

T o =0,l ° - 1 



with ^ JJJV (T C ) given by Eq. ©. Finally, Eqs. (l9ll"0l) 
become 

(6,J fc ;l/2,r 3 |Alg 1) ||a,J a ;l/2,r3) = LD^!±!l ( C 3) 
x E [AlxjA^+^idA0)]Ui\D^\jN) 

31 >3N 

(b, J b ; 1/2, T 3 \\M { e 0) \\a, J a ; 1/2, T 3 ) (C4) 
= \ E [^v(l)-^x^(0)](jilMsIiw) ■ 



With these modifications the possible consequences of the 
inclusion of more orbitals on the results of the different 
sections can now be considered. 

Section [IIII only concerns the form of the El operator, 
and does not depend on the assumed form of the wave- 
functions. 

Section IIV Al also is completely valid, as the consider- 
ations reported there do not depend on the details of the 
wavefunctions. 

Section IIV Bl depends on the assumed structure of the 
analogue and anti-analogue states. The choice given 
there presumably corresponds to an upper limit of the 
mixing. For example, in Eq. (|4ip. the choice of a pure 
.99 /2 orbit corresponds to the maximum possible value of 
the expectation value of r 2 /R 2 . Our conclusion, i.e. that 
this process is not able to explain the observed effect, is 
therefore even stronger if other orbitals are considered. 

It remains to consider Scction flV CI The sum on ji,jN 
must be included in Eqs. (|54U59p to obtain the one-body 
and the two-body contributions to the induced isoscalar 
El: 



1 



{b,J b ;\ 

(j N \(r/R) 3 \h) 



JN 



F (l)( 



J I 



3I,3N 



and 



(3N\{r/R)\ji) 



b ; J b ;l/2,T 3 \\M l °} b a ; J„,1/2,T E 



(C5) 



(C6) 



-C\ E ^iA^\{3N\rY^{f)\h) ((r/R) 2 ) v 



31, JN 



We obtain therefore 

_ (6, 7/2-; 1/2, T 3 \\M^\\a, 9/2+; 1/2, T 3 ) 

e ( J 3j = 



(C7) 



(6, 7/2-; 1/2, T 3 \M% 1] \\a, 9/2+; 1/2, T 3 ) 
= (_l)V2+r, 3C x 

gj, ,i„ (jw I ( T / R ) 3 \3l) hi J» A 3r ,3* ( X ) - A H Jn (0)] 

E jl<jN (3N\r/R\h) [A ]lM (l) + 3A lxU (0)} 
where f]jj,j N has the same meaning as in Eq. (|60l) . 
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